The evolution of GRB remnants 
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ABSTRACT 

The detection of the delayed emission in X-ray, optical and radio band, 
i.e. the afterglow of 7-ray bursts (GRBs), suggestes that the sources of GRBs 
are hkely to be at cosmological distances. Here we explore the interaction of 
a relativistic shell with a uniform interstellar medium (ISM) and obtain the 
exact solution of the evolution of 7-ray burst remnants, including the radiative 
losses. We show that in general the evolution of bulk Lorentz factor 7 satisfies 
7 oc t~°'* when 7^1, here at is mainly in the range 9/22 ~ 3/8, the latter 
corresponds to adiabatic expansion. So it is clear that adiabatic expansion is a 
good approximation even when radiative loss is considered. However, in fact, 
at is slightly larger than 3/8, which may have some effects on a detailed data 
analysis. Synchrotron-self-absorption is also calculated and it is demonstrated 
that the radio emission may become optically thin during the afterglow. Our 
solution can also apply to the nonrelativistic case (7 ~ 1), at that time the 
observed flux decrease more rapidly than that in the relativistic case. 

Subject headings: gamma-rays: bursts - radiation mechanisms: non-thermal 
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1. Introduction 

Since the discovery of 7-ray bursts nearly 30 years ago, their origin and emission 
mechanisms remain mysterious, a major reason is that their emissions other than X-ray 
and 7-ray have remained invisible. This situation has been changed dramatically since the 
launch of Italian-Dutch satellite BeppoSax. The delayed emission at X-ray, optical and 
radio wavelength has been detected from several GRBs due to an accurate determination of 
their position ( jUosta et al. 199"7a] , |Costa et al. 19975] ) The very long time afterglow (more 



than one month) strongly suggests that the sources of GRB are at cosmological distances. 

It is a natural prediction of the cosmological fireball model that 7-ray bursts should 
have an afterglow in X-ray, optical and radio bands ( [Meszdros fc Kees 1997] , [Faczynski fc 



Rhoads 1993 , Vietri 1997a ). After the main GRB event occurred, the fireball continues to 



propagate into the ISM, and thus relativistic electrons may be continuously accelerated 
to produce the delayed radiation on timescales of days to months. Several authors have 
discussed the radiation mechanisms of GRB afterglow and the agreement between the 
fireball deceleration model and the measurements is good ( [Waxman 1997a| , [Waxman 



1997b| , IWijers, Rees fc Meszdros 199?] ). However, these authors simply assume that the 



fireball expansion is adiabatic. On the other hand, Vietri (1997a; 1997b) has assumed that 
the fireball expansion is highly radiative and the surrounding matter is nonuniform, his 
results is also in agreement with the observations. 

In this Letter we explore the interaction of a relativistic fireball with a uniform ISM. In 
section 2 we derive the exact solution of the evolution of the bulk Lorentz factor 7 in general 
case, without assuming expansion adiabatic or highly radiative. Our solution can also 
extend to nonrelativistic case. In section 3 we discuss the emission features for adiabatic 
expansion and non-adiabatic expansion. The evolution of the synchrotron-self-absorption 
frequency is also calculated. Finally we summarize our results and give some implications 
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for future observations. 



2. The evolution of GRB remnant 

Whatever the sources of GRBs are, the observations suggest the fohowing scenario. 
Very large energy {E ~ 10^^ erg for cosmological distances) is suddenly released (T < 100 s) 
in a small region of space, so the initial energy density is so large that an opaque fireball 
forms ( Paczynski 1986 , [Goodman 1986 , Shemi fc Piran 1990| ), then the fireball expands 



outward relativistically. After an initial acceleration phase, the fireball energy is converted 
to proton kinetic energy. When the fireball is decelerated by the swept external matter, a 
strong shock can be created and electrons may be accelerated to very high energy. Then 
GRB is produced through synchrotron radiation or possible inverse- Compt on emission of 
these electrons. This occurs at decelerating radius ( [Meszdros, Rees fc Papathanassion 1994|) 



r. = i-^^^f" = 5 X 10^^(^)V3(^)-/3 cm (1) 
ri'^niiTLpC'^ rii 100 

where E = lO^^E^i is the total burst energy, rii is the density of ISM, rj = E/MbC^ 
represents the ratio of fireball energy to the initial rest mass energy of baryons, Mf, is the 
total mass of polluting baryons initially mixed with the fireball, and 70 is the Lorentz factor 
of blast wave at radius r^. It should be noted that at this time the shell thickness is similar 
to that of the heated ISM region, so the heated ISM and shell carry similar energy ( |Sari fc 



Piran 1995| , [Waxman 1997a|) . Therefore it is easy to show that 70 = ?7/2, and the ISM mass 



swept by the forward blast wave at this time is 2/// of the shell's mass, i.e. Mq/Mi, = 2/1]. 

Sari & Piran (1995) pointed out that the fireball energy dissipation may occur at two 
different places. They defined a radius te where the reverse shock becomes relativistic. If 
fd < fE, relativistic reverse shock can not be formed, and the shell loses its energy to the 
heated ISM at r^. Otherwise if > te, then most of the shell kinetic energy is converted 
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into thermal energy once the reverse shock crosses the shell (see also Waxman 1997a|) . 



This conclusion is valid only if the shell propagates with a constant width A. However, if 
the shell is expanding, then the shell width A = will not be a constant. According 
to the definition of Sari & Piran (1995), the parameter /(r) = n^^/rii = '-/P/rir^, where 
/ = ( — ^-^Y^'^ is the Sedov length. Then when the reverse shock becomes relativistic, i.e. 
l"^ I fijE) = 1, we can obtain 

TE = {^^f" (2) 

comparing with equation (1), we find that te ~ r^- So it is clear that when the fireball 
expanding, the density ratio / decrease with time. Initially 7^/ / -C 1 and the reverse shock 
is Newtonian. When r ~ r^; (r^) the fireball is decelerated by ISM and at the same time the 
reverse shock become mildly relativistic. 

After the GRB occurs, the fireball decelerates and a relativistic blast wave continues 
to propagate into the ISM to produce the relativistic electrons that may produce the 
delayed emission on time scales of days to months. The light curves of the afterglow 
are controlled by the deceleration of bulk Lorentz factor 7 and the intrinsic intensity 
In previous researches some authors ( | Waxman 1997a| , |Waxman 1997b| , |Wijers, Rees fc 



Meszdros 1997 ) assume that the fireball expansion is adiabatic, i.e. the radiative energy 



loss can be negligible, however some other people ( |Vietri 19974 [Vietri 1997b| ) assume that 



the expansion is highly radiative, i.e. the internal energy of the system is instantaneously 
radiated out, it is obviously that these are two extreme cases. In fact, the fireball should 
expand between these two cases, some fraction (not all) of the internal energy is radiated 
and lost from the system. So it is important to calculate the evolution of bulk Lorentz 
factor 7 when including radiative energy loss. We define to be the fraction of the internal 
energy emitted, then the total energy radiated per unit swept mass is given by ( [BlandforH 
fc Mckee 197"6| ) 



dE = -7(7 - l)esC^dM (3) 
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where M — ^nr^ninip is the swept ISM mass. For simphcity, we assume that the energy 
loss of the swept up matter is small, then the total kinetic energy can be written as 
E = (72 - l)Mc2. Thus we obtain 

dn_ ^ (1 + 6^7^ - e,7 - 1 

dM 2M7 ^ ^ 

which has the solution 

(7-1)^7 + 1/'^) ^ (.^ 
(70- 1)^(70 

where |/ = 1 + e^, 70 and Mq are the initial values of Lorentz factor and the mass of 
swept up matter at the radius Tq where the deceleration begins (ro — r^). When the 
expansion is highly relativistic (7 ^ 1), = i'M^)~^ — expansion is adiabatic 

(e^ = 0, y = 1), = (^)~^''^- Otherwise if expansion is extremely radiative (e^ = 1, y = 2), 
then ^ = (^)~^- However, it should be pointed out that, in the above calculation we have 
taken the value of as a constant, this is appropriate for describing the adiabatic and 
non-adiabatic limits (e^ = 0, 1), but for intermediate values the value of is almost certain 
not to be a constant, and in fact probably depends on gas parameters which depend on 
radius. 



3. the emission from GRB remnant 

The variation of bulk Lorentz factor 7 with observer time t can be obtained by 
combining relation dt — ^^^2^^^ and equation (5), where z is the cosmological redshift of 
GRB source. For 7 3> 1 (/5 ~ 1) it is given by 

- = (6) 
7o to 

where to = {1 + z)ro/2{l + 3y)7oC is the typical duration of GRB event, at = 3y/2(l + 3y). 
The fraction of the energy radiated can be estimated as = CeCsyn, where Cg represents 
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the fraction of the energy that is occupied by electrons, the typical Lorentz factor 
of electrons in the comoving frame can usually be expressed as 7e = ^e^il — l)j so 
Cg — = 1 (ce = I) represents the energy equipartition between electrons and 

protons. The synchrotron radiation efficiency can be expressed as egyn — .-i"""-! (we 
assume that synchrotron radiation is the main mechanism for GRB afterglow), where 
tsyn = QT^trieC/aT'yeB^ ^nd = r/^f3c are the synchrotron cooling time and expansion 
time in comoving frame respectively, B — (^^8717(7 — l)nimpC^)^/^ is the comoving frame 
magnetic field, me{mp) is the mass of electron (proton), ut is Thomson cross section. 
Therefore we expect should lie in the range ~ 1/2, and y in the range 1 ~ 3/2, so at 
should between 3/8 (adiabatic expansion) to 9/22. Thus it is somewhat unexpected that at 
lies in a very narrow range, for general case at is only slightly larger than 3/8 (adiabatic 
case) , so it is reasonable to think that adiabatic expansion is a good approximation and the 
radiative energy loss could be treated as a small correction. Therefore in the following we 
first consider the emission features at adiabatic expansion. 

3.1 adiabatic case From equation (5) we see that for adiabatic case {y = 1) the bulk 
Lorentz factor 7 can be expressed as 

7 = (l + (7o^-l)^)^/^ (7) 

without any assumption. Here Tq is the initial value of the deceleration phase, so Tq = r^. 
Prom equation (1) = 2E /rfninipC^ = E / = /^/27q, so ^y^r^ = 1^/2. Thus for 
7o ^ 1 we have 

where x = r/fc, fc = //2^/^, for 7 3> 1 from equation (6) we obtain 

7 = 278(1 + zf'Ell'n^'/'r'/' (9) 

This expression is valid only for relativistic expansion, it breaks down for 7 ~ 1. The time 
through which the blast wave remains relativistic is about a month. 
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If synchrotron radiation is the main mechanism to produce the GRB afterglow, we can 
estimate the time delay between the beginning of the afterglow and the onset of the optical 
flash. The observed photon energy of synchrotron emission from typical electrons is (7 ^ 1) 

= = ' + ^) 

where = 4.413 x 10^^ G is the critical magnetic field. Then combining this equation with 
eq.(9), we have 

^ 1.1 (1 + zy/'Eli^^r'/'^ej^el' (days) (11) 

Thus considering the uncerntain of the parameters, the value of tm (optical flash) could 
range from a few hours to about 3 days. 

It can easily be seen from equation (10) that the typical photon energy cannot possibly 
enter the radio region when 7 ^ 1, so we expect that until the fireball expansion is medium 
relativistic does the synchrotron emission peak at radio band. It is interesting to note 
that our solution (eq.(8)) can be extended to non-relativistic case. From eq.(8) the shell 



velocity (3 = l/vT+x^, then the relation between observer time and radial distance is 
t = f ^^+^^{3^"^)'^" = JiVTT^ - l)dx, so when a; < 1 (7 > 1), r oc while when 

X ^ 1 (7 ~ 1), r oc 

The light curve of afterglow is also dependent on the intrinsic photon 
spectrum. We assume the comoving intensity is oc for e < and 

oc for e > e^. From equation (8), we can see, in the relativistic case, 
7e oc i? oc 7, Cm oc 7e-B7 OC t"^^"^,!'^^ OC UeBAr oc t^^^^, so the observed 

peak flux F^^^^ oc t'^'j^I'^^ oc ~ constant ( |Meszdros fc Rees 1997| ) Before crosses 
the optical (radio) band, the flux oc -F£„(^)" oc t^"/^, while after crosses the 
optical (radio) band, we have F^ oc F^^{-^)^ oc t^^^"^. In the non-relativistic case, 
7e oc t~^^^, B oc t"^/^, e-m oc t"^, /^^ oc and F^^ oc t^^^. So the observed optical 

(radio) flux for e < is oc F.(^r oc for e > is oc F.(^y oc 
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These features have also been discussed by Wijers, Rees & Meszdros (1997). Fig.l 
give an example of the variation of observed flux with time t. We have taken 
a = 0, jS = —1, E = 5 X 10^^ ergs. The solid, dashed and dot-dashed lines correspond to 
optical, x-ray and radio flux respectively. It can be seen from Fig.l that the slope (e > e^) 
changes from -1.5 (relativistic case) to -2.4 (non-relativistic case) gradually. 

3.2 non-adiabatic case It can be seen from equation (5) that for non-adiabatic case 
(?/ > 1) we cannot obtain the explicit expression for the evolution of bulk Lorentz factor 
7, only when 7 ^ 1 it can be written as the form of eq.(6). In order to compare with 
the adiabatic case, we write at = 3/8 + a', where a' should be less than 3/88. Then we 
can write the evolution of Lorentz factor as 7 = 7ad(^)~"', (~ 1 ~ 100 s) is the typical 
duration of 7 ray flash, then we see that after one day 7/7ad ~ 0.75 for to ~ 10 s, a' = 3/88, 
the exact value will depend on the parameter e^. In the meanwhile, from eq.(ll) the time 
delay between the GRB event and the onset of the optical flash will be shorter by a factor 
(t/to)~^"' ~ 0.43, i.e. it could be range from a few hours to more than one day. 

This will also affect the light curve of afterglow. Just as discussed above, for relativistic 
case, 7e cx oc 7, and oct (2+^"'). It should be noted that, in the non-adiabatic 

case, the electron cooling time is shorter than the expansion time, and the effective width of 
the emission region is narrow than the width of swept up matter, so /^^ oc UeBctsyn oc ti"*""' 
( [Meszdros, Rees fc Wijers 1997] ), and the observed peak flux F^^ oc t^~'^°'' is not a consant, it 
increases with time. Then the flux oc tih+^°')+Mo'-i)a' t < e^n and F^ oc 
for e > e^. Therefore we see that for t < t„i the flux F^ may either increase or decrease 
with time, which depends on the values of a and a' . 

For non-relativistic case the variation of 7 with t is more complicated. However, it 
should be pointed out that radiative energy loss may be important only in early time. The 
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ratio of synchrotron cooling time to the expansion time in the comoving frame is 

^ = 7 X lOVCTsWto'ih''''-' (12) 

lex I-O 

Taking at = 9/22, for Uyn/te. = 1 we obtain ~ OAee'^'^B^' Eli'nf day. 
Therefore we expect that the synchrotron coohng time to be longer than the dynamical time 
for long delay afterglow (usually te < 1 day). When tgyn > tex, the synchrotron radiation 
efficiency egyn decreases rapidly and so -C 1, at that time the adiabatic expansion is a 
very good approximation. 

3.3 synchrotron-self-absorption We have shown that the typical photon energy 
em might be impossible to enter the radio region when 7 ^ 1, so the detected radio 
flash within a few days after GRB ( [brail et al. 1997 ) may be possible due to the source 



becoming optically thin. As Waxman (1997b) and Vietri (1997b) have pointed out, 
the synchrotron-self-absorption optical depth for photon energy e < may scale as 
T = Tm(^)~^ due to the presence of a low energy electron population. The optical depth 
Tjn can be estimated as ~ 1.2 x 10"-'^^n|;^^^^^^^^^^7^/^(7 — l)^^^^'^(3t, then for adiabatic 
approximation, the absorption frequency for 7 ^ 1 is Uab — "^fh^"^ ^"^^"^ EK^f^^y^ GHz, 
while for 7 ~ 1, — 0.2nl^^^~^^'^^^H]^^y GHz. So the absorption frequency first decreases 
with time as Uab oc while when the expansion is non-relativistic, Uab oc t^^^. When 

photon energy e > e^, t oc e~^~^^^^'^, where p is the index of electron distribution, then 
i^ab = Taking p = 3, then u^b = 317nfej'^T'EU'%^^^'^ GHz for 7 » 1, and 

Uab = 1556nl^^'^^j^'^^^^'^El(^t'^^yGHz for 7 ~ 1. The observed fiux may consists of several 
components. If Uab < i^m, then F, = F,„(^)"(-^)-+2 for u < u,b, F, = F,^{^Y for 
yab<y< Vm. and F, = F,^{^Y for v > u^. If i^ab > then F, = F,„(^)/3(-^)5/2 fo^. 
u < Uab, and F^ = -Fy,„(— )^ for u > Vab- These components can be seen in our calculated 
radio fiux in Fig.l. 



4. Discussion and conclusion 



The detection of 7-ray burst in the optical and radio bands has greatly furthered our 
understanding of the objects. Some people have discussed the afterglow emission using the 
adiabatic expansion or highly radiative expansion. Here we have investigated the evolution 
of the fireball Lorentz factor 7 when it interacts with the ISM. We find that even though 
including the radiative energy loss, the difference of the evolution of the 7 between our 
results and that in adiabatic case is small. So adiabatic expansion can be treated as a good 
approximation. 

From eq.(6) wc can also obtain the evolution of the total system energy (when 7 ^ 1) 
E/Eq — {t/tQ)~^^''/^'^^^^''\ which is different from the results of Sari (1997), who obtained 
E (X ^-I'^'e^/ie. This is because Sari had used the self-similar solution to obtain his results, 
while our results are obtained by solving the energy loss equation without any assumption. 
Our results show that the decrease of total energy is slower than that predicted by Sari 
(1997). 

We have shown that the energy loss due to radiation may be important only in early 
time of afterglow (usually t < 1 day), since the synchrotron cooling time may be much 
longer than the expansion time for long delay emission, so we argue that the detection of 
afterglow within short time after the GRB event is very important, since it may provide 
more information about the GRB sources, the radiation processes, and the surrounding 
matter features, etc.. 

The afterglow has been detected for more than one month, the interesting question is 
how long can the afterglow be sustained? Prom Fig.l we see that the optical flux decays 
with time as a power law F oc with the index n increases gradually from |/? to (| + 3/3). 
Therefore we argue that the optical flux should be detected for longer time without suddenly 
cutoff. In addition it is also possible to observe radio emission for very long time. 
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The detection of very long afterglow strongly suggests that the sources of GRBs are 
at cosmological distances, which mean the total energy of GRB event should be about 
]^g5i-52 Qj.gg_ Now the most fundamental problem, the ultimate energy source and the 
physical processes leading to the fireball formation has not yet been solved. Recently the 
popular model-coalescing of two neutron stars- seems to be difficult to account for the 
GRBs ( |Ruhert et al. 199 7] ), and now a new way, very strong magnetic field {B ~ 10^^ G) 
combined with rotation, has been suggested to power the fireball ( [Paczynski 19971 ). We 
expect the detection of GRB afterglow can provide information about the origin of GRBs. 
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Fig. 1. — The evolution of observed flux with time t. We have taken a — 0, (3 — —1, E — 
5 X 10^^ ergs. The sohd, dashed and dot-dashed hues correspond to optical, X-ray and radio 
flux respectively. 



